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Introduction 

The awareness that transmission dynamics of infectious diseases, the preventive and 
treatment strategies that are administered follow certain laws that could be modeled 
mathematically have been a beautiful child to the human conscience since the first incidence 
of epidemic outbreak [1]. 

The awakened consciousness that preventive and treatment strategies which are 
administered in the approach to finding solution to an epidemic outbreak require some 
mathematical calculations have contributed immensely in removing the subterfuge mist 
that has surrounded the elucidating of the transmission dynamics of infectious disease, the 
control strategies and the prevention of epidemic from invading a population of susceptible 
and has led us out into the light of strategic epidemiological control measures [2]. 

The mathematical computations which have been involved in the study of the dynamics 
of infectious diseases have brought a great surge of excitement in the hearts of modelers, 
governments, public health workers and all stakeholders as it has created awareness with 
regards to transmission dynamics and control measures of infectious diseases [3]. 

Concerning the malaria disease, (WHO, 2016) reported a global estimate of 212 million 
cases, of which the most occurred in Africa, 429000 malaria deaths, of which 303000 
occurred in children under 5 years of age. These estimates rank malaria to be virulent 
disease and one of the top two killers among vector borne diseases. The disease is caused 
by parasites 

that are transmitted through the bites of a host mosquito. An infected person with the 
disease may experience loss of appetite, headache, fever and fatigue. The disease has an 
adverse effect on the well-being of infected person and the task of achieving 

malaria-free world throws a challenge to the public health organization which 
requires a strenuous and strategic efforts. The burden of the disease includes; fertility loss, 
premature mortality and loss of savings through treatment and medical cost, loss in worker 
productivity due to weakness in body. Malaria recede growth in population by means of 
morbidity and mortality in especially children. The disease poses threat to the development 
of fetal in the early stages of pregnancy in women as a result of loss of immunity. Control 
measures of the disease include the preventive intervention of sleeping under insecticide 
treated nets and spraying the breeding environments of the vector by insecticide spray and 
the treatment measure of seeking medical attention in the early stage of the disease [4]. 

Notwithstanding, there is an urgent need for better understanding of the disease 
transmission and development of effective optimal therapeutic strategies for the prevention 
and control of the disease due to drug and insecticide resistance challenges posed by the 
parasite [5]. 
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Mathematical modeling has been a fruitful ground out of which 
have grown control therapeutic strategies for the malaria control into 
fruition. The contributions of mathematical modeling to public health 
organizations in their intervention strategic measures of controlling 
malaria disease have been overwhelming [6]. Abdullahi et al [7], 
proposed a deterministic malaria model that explains the dynamics 
of the disease and measures the effectiveness of anti malaria drugs 
in the infected compartment. Mandal et al [8], presented a review 
of existing mathematical models on malaria disease and their 
contributions with regards to the explanations of the transmission 
dynamics. Their work elaborates the evolution of modeling strategies 
to describe malaria incidence with the inclusion of critical features 
of host-vector-parasite interactions, by using an already existing 
mathematical model. Pariyaporn et al [9], used two malaria models: 
malaria model with standard incidence rate and malaria model with 
nonlinear incidence rate to explain the effect of choice of incidence 
function that results in occurrence of backwards bifurcation in these 
models. Feng et al [10], presented a deterministic model of malaria 
transmission with standard incidence rate and treatment, which was 
used to numerically stimulate the human data of the disease reported 
by the Chinese ministry of health from 2002 to 2013. 

Besides, optimal controls models have come up with astonishing 
models coupled with various explanations of the transmission 
dynamics of the vector-parasite-host and have induced control 
therapeutic preventive strategies which could be administered 
for the prevention of the disease as well as treatment strategies 
for the treatment of the disease at a minimal cost [11, 12, 13]. Kim 
et al [14], presented a Plasmodium Vivax malaria transmission 
model that includes control terms by using a deterministic system 
of differential equations to analyze the model mathematically and 
numerically. Okosun et al [15], examined the conditions under which 
it is optimal to eradicate malaria disease by proposing and analyzing 
a deterministic model for the transmission of malaria disease that 
includes mass action. Fatmawati and Tasman [16], formulated an 
optimal control model of malaria transmission that considers the 
resistance of malaria parasite to the anti- malaria drugs. The model 
examined the impact of vaccination and treatment strategies on the 
disease transmission. Okosun et al [17], investigated the effectiveness 
and cost-effectiveness of malaria model with three preventive 
measures: the use of treated bednets, spraying of insecticide and a 
possible treatment of infective humans that blocks transmission to 
mosquitoes. Okosun and Makinde [18], examined the synergistic 
relationship between malaria and cholera disease in the presence 
of treatment by formulating a co-infection mathematical model of 
malaria-cholera. Mwanga et al [19], studied the optimal control 
practices of malaria. The model addresses the implementation of a 
catalog of optimal control strategies in the presence of uncertainties 
of parameter, which is generally a typical problem of infectious 
disease data. Rafikov et al [20], investigated the development of 
transgenic mosquitoes which may provide a new challenge and 
serve as an effective weapon of diseases control as a result of their 
resistant to diseases. Lashari et al [21], assessed the impact of some 
antimalaria control measures by reformulating an existing model 
as an optimal control problem. Agusto et al [22], applied optimal 
control theory to a deterministic system of differential equations to 
investigate an effective optimal strategy for controlling the spread 
of malaria disease using treatment, insecticide treated bednets and 
spray of mosquito insecticide as the system control variables. 

Enlisted in these models are useful and successful companions 
and comrades in elucidating the transmission dynamics of the disease, 
deducing of control intervention measures to be embarked on by all 
stakeholders as well as insinuating imperative therapeutic strategies 
for medical decision making, yet there has not been a model that has 
been able to answer all clinical and biological information that are 
empirically observed regarding the disease as well as incorporating 
all important factors and variables in the model due to the complexity 
of the disease. In this research article, we consider the malaria 
model with standard incidence rate as proposed by [9], as studied 
in section 2. The section 3 deduces an optimal control problem that 


assesses the impact of some control strategies: prevention, treatment 
and application of insecticide spray on the vector by incorporating 
time dependent control functions. The necessary conditions for an 
optimal and the corresponding states are then derived by employing 
the Pontryagin's Maximum Principle. Finally, in section 4, the 
resulting optimality system is numerically solved and computed to 
investigate the optimum control strategy that would be efficacious 
to be implemented in reducing the number of exposed and infected 
humans. Further, in section 6-7, Stochastic version of system (1) is 
obtained by introducing a random perturbation to the equations 
of the systems. In section 8, stochastic numerical method-Euler- 
Maruyama is discussed. Section 9 gives numerical examples to show 
the dynamics of the deterministic, control and the stochastic systems, 
and to illustrate the differences in the dynamics of the models. 

Model Framework 


This section presents a malaria model with standard incidence 
rate. The human population at time t is categorized into four 
classes: susceptible human S H (t), exposed human E H (t), infectious 
human I H (t), recovered human R H (t), respectively. The mosquito 
population is divided into three classes: susceptible mosquito S y (t), 
exposed mosquito E y (t) and infectious mosquito I y (t), respectively. 
The total populations of human and mosquito at any time are given 
by N H (t) = S H (t) + E„(t)+ I H (t)+ R H (t) and N v [t)= S v (t)+ EJt)+ I v (t) 
respectively. The mosquito population does not include immune class 
as mosquitoes never recover from infection, that is, their infective 
period ends with their death due to their short life-cycle. Further, it 
is accepted that the recovered human have some level of immunity to 
the disease and do not get clinically ill, but still harbor low levels of 
parasite in their blood streams which are transmitted to a susceptible 
mosquitoes during bites [9]. The mathematical differential equations 
of the dynamics of the malaria model as proposed by Pariyaporn et 
al [9] are: 
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The model assumed that people enter the susceptible class either 
through birth or immigration at a recruitment rate^-A- When an 
infectious mosquito bites a susceptible human, there is some finite 
probability that the malaria parasite will be passed on to the human 
at the rate of o v p H ^-, and the person will move to the exposed class. 
Individuals from the exposed class enter the infectious class at a rate, 
that is the reciprocal of the duration of the latent period. Recovered 
infectious humans move to the recovery class at a rate and die from 
the infection at a rate. It is assumed that recovered individuals have 
temporary immunity that can be lost and are again susceptible to 
re infection at a rat ePH. All human classes leave the population 
through the same natural death rate . The infectious human 
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leaves the population at a per capita disease-induced death rate. For 
the mosquito population, susceptible mosquitoes are recruited by 
birth at a constant rate, independent of the actual number of adult 
mosquitoes. This assumption reflects reality, since only a fraction of a 
large reservoir of eggs and larvae matures to the adult stage, and this 
process does not depend directly on the size of the adult mosquito 
population. All mosquito populations die due to their finite life span 
at the same rate of. Susceptible mosquitoes are infected by biting 
the infectious human and recovered human at the infection rate 

°v f ^ vI " + J 3v,,R " I Susceptible mosquitoes which are infected move to 


the exposed mosquito class and progress to the infectious mosquito 
class at a rate v v . 

Optimal Control Strategies 

In this section, the state system (1) is modified to assess the 
impact of some control strategies: prevention, treatment and 
application of insecticide spray on vector. In the human population, 
the associated force of infection is reduced by a factor of (1 - CO x (t) 
) where CO x (t) represents the use of preventive control of sleeping 
under treated mosquito net to prevent direct contact and bite from 
infected mosquito. ^ ( t ) re P resents an application of insecticide 
spray on the bleeding grounds of the vector. It is undeniable fact that 
mosquito bleeds and flourish under a favorable climate conditions 
such as rainy seasons and hot climatic conditions [23]. Therefore, 
strenuous effort is required in application of insecticide at the 
bleeding fields of the vector at these times. Thus, the reproduction 
rate of the mosquito population is reduced by a factor of (1 - ^ (7) 
). The control function C0 3 (7) measures the rate at which infected 
individuals are treated at each time with malaria drugs for the 
control of the disease. We assume that ^ / individuals at any time 
(t) are removed from the infective class 3 ancl added to the removed 
class. Further, we assume that the mortality rate of the vector 
population increases at a rate proportional 8 X I y C0 3 {t ), where § 

> 0 is a constant rate. With regards to these underlying assumptions, 
an optimal control model for malaria disease is formulated that 
deduces prevention and treatment strategies with a minimal cost of 
implementation. Hence, the dynamics of system (1) are modified to 
the following system of equations: 
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The objective of our work is to minimize the number of exposed 
and infected human individuals and the infected vector population 
and maximizes the total number of recovered individuals through 
preventive and treatment strategies, by employing feasible minimal 
time dependent control variables CO x ( t ), C0 2 ( t) and C0 3 (t) 
respectively. 


With appropriate initial conditions, we consider an optimal 
control problem to minimize the objective functional given by 


J(w 1 ,w 2 ,w 3 ) = \^A t E H +A 2 I h + A 3 I v + a 2 w 2 2 + a 2 w 2 )j dt (3) 


The quantities A v A 2 and A 3 denote the weight constant of 
the exposed and infected human individuals and the infected 
vector population. Again, the quantities w v w 2 and w 3 are weight 
constants for minimizing the number of exposed and infected human 
individuals and treatment of infected human individuals. Further, the 

terms 1 „ 2, 1 „ 2 and ^ n 2 represent the cost associated with 

-d x Wi -a 2 w 2 -a 3 w 3 ^ 

the minimizing the exposed and infected and treatment of infected 
human individuals. 

We choose a quadratic cost on the controls as a reflection of what 
is in other literature on epidemic control models [24]. Here, we seek 

an optimal control W l , W 2 , W 3 such that 

J(m\,w* 2 ,w* 3 ) = min{7(w 1 ,w 2 ,w 3 | w 1 ,w 2 ,w 3 gW)} (4) 

Where 

{w l ,w 2 ,w 3 1 w f (Ois lebesgue measurable with, 0 < w f (0 < 1, i=l,2,3} (5) 


Applying the Pontryagin's Maximum Principle [25], the system 
(2) and (3) are converted into minimizing the Hamiltonian H, with 
respect to w v w 2 and w 3 
Where 


H (S H , E h , I H , R h , Sy , Ey , Iy , W, , Aj, A 4 , A 5 , 2 6 , Ay, t ) — A^E H + A^I H + A^I V 
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Where /t ; for i=l,.,7 are adjoint variables to be determined 

* * * 

Theorem 3.2: There exists an optimal control W* = ( W l , W 3 , W 3 
)EW such that 


J(w x , w 2? w 3 ) /(w,,w 2 ,w 3 ) 

subject to the control system (2) with initial conditions at t=0 . 
Proof: The existence of an optimal control is proved by employing 
the result in [26]. We note that the control and state variables 
are nonnegative values. Hence, the necessary convexity of the 
objective functional in w v w 2 and w 3 are satisfied in this minimizing 
problem. The set of all control variables ( w i ? w 2 ? w 3 ) Gw is also 
convex and closed by definition. The optimal system is bounded 
which determines the compactness needed for the existence of the 
optimal control. In addition, the integrand in the functional (3), 

A x E h + A 2 I h + A 3 I v +— +a 2 w 2 +a 3 w 3 ) j s CO nvex on the control 

set V. Also we can easily see that, there exist a constant V>1 and 
positive numbers V v V 2 such that 

J( Wl, W 2 , w 3 ) > V l (| w, | :2 +1 w 2 \ 2 +1 w 3 1 2 ) 2 - V 2 


because, the state variables are bounded, which completes the 
existence of an optimal control. In order to find the optimal solution, 
we apply Pontryagin's Maximum Principle [27] as follows: 

Given that (x,w) is an optimal solution of an optimal control 
problem, then there exists a non trivial vector function x=(x i , x 2 , 

x 3 .x ) which satisfies the inequalities 

dx _ dH (t, x, w, X) 
dt dX 

dH(t,x,w,, l) (7) 
dw 

^ _ dH (t, x, w, X) 
dx 

Now we apply the necessary conditions to the Hamiltonian H in (6). 
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Theorem 3.3: Given that g* f 9 S* ,E*, I*) are optimal 

( * * *. 

W l , W 2 , W 3 ) are associated optimal control 
variable for the optimal control problem (2)-(3), then , there exists 
adjoint variables xi, for i=l,2,.7, which satisfies 

, ( Err + I it + Rtt ) ( J 6 I it + P v fjRfJ ) 
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with transversality conditions (or boundary conditions), 
A .(^)=0,i=l,2,.,7 (9) 

* * * 

Furthermore, optimal control 14^ , W 3 , 1V 3 are given by 
w; =max jmin ,l|,o| ( 10 ) 

a t N H ’ J’ J (11) 

w* = max jmin 7 ^/* + (^ - /l 4 )/^, 'j’ 0 } (12) 


w 7 = max < mm 


Proof: The adjoint equations and the transversality conditions are 
determined by employing the Hamiltonian (6). By putting S H = S* H 
(t), E h = e‘ h (0 ,1 = r H (t), R = R‘ h (t), S„= s; (t), E„= e; (t) and 

I y = I* v and differentiating the Hamiltonian with respect to S Rt E H , I H , 
R h , S y , E y and I v respectively, we obtain (8). By solving the equations 
dH 

M =0> dH = q, ^L = Qon the interior of the control set and 

dw x dw 2 dw 3 

using the optimality conditions and the property of the control space 

W , we obtain (10)-(12). Here, we empty the formulas (10)-(12) for 

* * * 

w*= ( W l , W 2 , FV 3 ) the characteristic of the optimal control. The 
optimal control and the state are found by solving the optimality 
system, which consists of the state system (2), the adjoint (8), initial 
conditions at t=0, boundary conditions (9), and the characterization 
of the optimal control (10)-(12). To solve the optimality system, 
we use the initial and transversality conditions together with the 
characterization of the optimal control (14^ , 4V 2 , W 3 ) given by 
(10)-(12). In addition, the second derivative of the Langragian with 
respect to wl, w2 and w3, respectively, are positive , which shows 

that the optimal problem is minimum at controls W ] , W ? and VR, . 

* * * 

By substituting the values of 14^ , W 2 and W 3 into the control 
system, we deduce the following system of equations. 
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dt dt dt dt dt dt dt 

Solving numerically the above systems (13) and (14) gives the 
optimal control and the state. 


Numerical examples and Discussion 

In this section, we asses by means of numerical approach 
the effect of the control strategies on the transmission dynamics 
of malaria model. In order to achieve this, an iterative scheme is 
applied in solving the optimality system: state system and adjoint 
system. First, the state system of equations is solved with arbitrary 
guess for the controls over a simulated time frame by employing 
fourth order Runge-kutta scheme. As a result of the presence of 
boundary conditions, the adjoint system is solved by backwards 
fourth order Runge-kutta using the immediate iterative solutions of 
the state equation. The controls are updated by means of a convex 
combination of the previous controls as well as the characterizations 
(10), (11) and (12). The entire process is repeated until the values of 
the unknowns at the previous itearations are closed to the one at the 
current iteration. 


The models assess the impact of control strategies on the 
transmission dynamics of Malaria disease with standard incidence 
rate. We investigate the effect of control strategies: prevention, 
treatment and application of insecticide spray on vector. This is done 
by numerically simulating the results of the dynamics of the disease 
to ascertain the stated scenarios of the strategies with parameter 
values S H (0) = 950000, e h (0) = 500000, 4 (0) = 25000, R H (0) = 
10, S v (0) = 850000, E v (0) = 300000, I v (0) = 850000, a t = 500, a 2 
= 30, a 3 = 10, A t = 10, A 2 = 10 and A 3 = 10 from [9], 

We further assumed that the weight factor, a 1# associated with 
control w 1 is greater than A v A 2 , A 3 , a 2 and a 3 respectively, which are 
association of controls w 2 and w . This holds on the ground that the 
cost of implementing w 1 includes, educational campaign programs in 
our public institutions such as schools, health centres and churches 
and media advertisements such as Television and radio educational 
campaigns to educate the public on the use of mosquito nets and the 
danger of exposing oneself to the mercy of the vector especially in 
pregnant women and infants. The cost of applying insecticide spray 
includes the insecticide cost and labour cost while treatment cost 
includes hospitalization, medical examination and the administration 
of malaria drugs in the infected persons. 

The parameter values used in the simulations are estimated based 
on the malaria disease as provided in Table 1. Other parameters were 
chosen arbitrary for the purpose of the numerical simulation. 

Figures 1-2 represent the number of susceptible Human [S H ) 
without and with controls for a 1 = 500. In the absence of control, the 
susceptible (solid curve) decreases gradually in the first 30 days and 
then rises slowly until almost the susceptible population is preserved. 
This is due to the fact that in the absence of control, more susceptible 
become infected with the disease and lose their status of being 
susceptible, until medical treatment has been sought. In the presence 
of controls, the susceptible (dashed curve) decreases sharply in the 
first 15 days and then experience a swift rise in the population, until 
almost their population are maintained in about lOOdays of the 
infection. This is possible due to the negligence, errors and difficulty 
of implementing these control measures, thereby exposing more 
susceptible to the vector, causing a lot of casualties in the early stages 
of the pandemic. But when these measures are properly instituted 
and monitored, the situation is reversed, thereby experiencing a swift 
rise in the curve. 
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Parameter 

Description 

Estimated value 

Reference 

n H 

Recruitment rate of human 

2.5 

[9] 

Ph 

Transmission probability from an infectious 
mosquito to a susceptible human 

0.9 

[9] 

A 

Transmission probability from an infectious 
human to a susceptible mosquito 

0.8 

[9] 

PvH 

Transmission probability from a recovered 
human to a susceptible mosquito 

0.009 

[9] 

Ph 

Natural death rate of human 

0.00004 

[9 

S„ 

Death rate from the disease 

0.00354 

[9] 

Yh 

Infectious human recovery rate 

0.003704 

[9] 

v H 

Progression rate from E H to I H class 

0.08333 

[9] 

K 

Progression rate from E v to I v class 

0.1 

[9] 

Ph 

Rate of loss immunity in human 

0.0146 

[9] 


Recruitment rate of mosquitoes 

500 

[9] 

<J V 

Biting rate of vector 

2.9 

[9] 


Natural death rate of vector 

0.0714 

[9] 

a h , ol v 

The level at which the force of infection saturates 

(0,1) 

[9] 


Table 1: Description of variables and parameters of the Malaria Model (1) 
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Figure 2: The plot represents population of susceptible Human 

with control 


Similarly, figures 3-4 represent the number of Susceptible 
Mosquito ( S v ) without and with controls for a 1 = 500. When there 
are no controls, the Susceptible Mosquito (solid curve) decreases 
sharply in the first 50 days until all the susceptible population are 
infected with the disease and leaves no population of susceptible. This 


is due to the fact that in the absence of any effective control measure, 
the mosquito-human contact is not prevented and the vector gets 
contact with human easily and feeds on his blood, thereby causing 
more people to be infected with the disease. However, in the presence 
of controls, the mosquito-human contact is minimized to some extent 
causing fewer people to be infected with the disease. Hence, there is 
a gradual decrease in the control susceptible mosquito curve (dashed 
curve) for the first 70 days until all the population gets infected with 
the disease. 

Figures 5-6 represent the number of Exposed Human ( E H ) 
without and with controls for a 1 = 500. Evidently, in the absence of 
control, the mosquito-human contact is not prevented. Hence, more 
people are expected to be exposed to the disease. This is shown in 
the sudden rise of the exposed curve without control (solid curve) in 
the early stages of 20 days and rises until almost all the population 
becomes exposed to the malaria disease. In the presence of controls, 
the mosquito-human contact is minimized drastically and people 
exposed to the disease eventually gain immunity through treatment 
or stronger immune system. Hence, the curve with control (dashed 
curve) experiences a gradual decrease for the first 60 days and 
decreases until all the exposed human leaves the population. 

Further, figures 7-8 represent the number of Exposed Mosquito 
( E v ) without and with controls for a 1 = 500. In the absence of 
control, the exposed mosquito (solid curve) decreases gradually for 
the first 50 days until all the exposed population is infected and leaves 
no population of exposed. In the presence of controls, the exposed 
mosquito (dashed curve) decreases sharply, and their population 
are maintained until about 30 days where all their population 
degenerated due to being infected. This is due to the fact that in the 
presence of controls, the mosquito- human contact is minimized and 
this prevent more mosquitoes from getting exposed, and enhance the 
exposed mosquitoes from leaving the population at a faster rate. 

Figures 9-10 represent the number of Infected Human ( I H ) 
without and with controls for a 1 = 500. It can be ascertained that, in 
the absence of controls, the mosquito-human contact is not prevented. 
More humans would be infected when they are bitten by mosquitoes. 
This is reflected in the infected human curve (solid curve). In the 
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Figure 3: The plot represents population of susceptible Mosquito 
without control 
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Figure 4: The plot represents population of susceptible Mosquito 
with control 






Figure 6: The plot represents population of Exposed Human with 
control 



Figure 7: The plot represents population of Exposed Mosquito 
without control 
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Figure 8: The plot represents population of Exposed Mosquito 
with control 



Figure 9: The plot represents population of Infected Human 
without control 
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Figure 10: The plot represents population of Infected Human 
with control 
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absence of controls, the curve rises sharply in the first 10 days and 
rises until it maintains its level for the rest of the days, until almost 
all the population becomes infected. The presence of controls is 
witnessed by an unprecedented rise in the infected human curve 
(dashed curve) which however reverse within 5 days and decreases 
sharply for the next 55 days where the infected human degenerates 
and leaves no population of infected. This suggests that the control 
strategies proposed were effective in minimizing the infected human 
population drastically. 

Similarly, figures 11-12 represent the number of Infected 
Mosquito ( Iy ) without and with controls for a 1 = 500. In the absence 
of controls, the infected mosquito curve (solid curve) decreases 
slowly and the infected mosquitoes are not completely removed from 
the population at the end of 100 days. 

However, with the introduction of controls, the infected mosquito 
curve (dashed curve) decreases sharply and the infected mosquitoes 
are wiped out from the population before 70 days. This insinuates 
that the control strategies proposed have been effective in minimizing 
the infected mosquito population drastically. 

Figures 13-15, present an optimal control effort to educate the 
public of the use and importance of implementing the preventive 
control of sleeping under treated mosquito net to prevent direct 
contact and bite from infected mosquito, w v the application of 
insecticide spray on the bleeding grounds of the vector, w 2 and the 
rate at which infected human seeks medical treatment, w 3 . The plots 
indicates that the preventive control w was at the upper bound 
for = 72, before it decreases gradually to the lower bound at = 100, 
while the optimal preventive control w 2 and treatment control 3 both 
maintain the maximum of 100% for = 99 before they drop sharply to 
the lower bound at = 100. This implies that minimum effort would 
be required for implementing the application of the insecticide spray 
on the bleeding grounds of the vector and treatment of the infected 
human than implementation of the preventive control of sleeping 
under mosquito net at a 1 = 10000. 




Figure 12: The plot represents population of Infected Mosquito 
with control 
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Figure 13: The plot represents Optimal control w 1 with <^=10000 
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Figure 14: The plot represents Optimal control w 2 with 0^=10000 
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Figure 15: The plot represents Optimal control w 3 with 0^=10000 


Similarly, figures 16-18 represent the control effort of 
implementing the preventive control of sleeping under treated 
mosquito net, w 1# the application of insecticide spray, w 2 and the 
rate at which infected human seeks medical treatment, w 3 . Here, 
we observed that the preventive control w was at the upper bound 
for = 98, dropping gradually to the lower bound at = 100, while the 
optimal preventive control w 2 and treatment control 3 both maintain 
the maximum of 100% for = 99 before they drop sharply to the lower 
bound at = 100. 

This suggests that minimal effort is required to implement the 
application of the insecticide spray on the bleeding grounds of the 
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vector and treatment of the infected human than implementation of 
the preventive control of sleeping under mosquito net at a 2 = 10000. 

Finally, figures 19-21 represent the control effort of implementing 
the preventive control of sleeping under treated mosquito net, v\q, 
the application of insecticide spray, w 2 and the rate at which infected 
human seeks medical treatment, w 3 . Again, we see that the preventive 
control w x was at the upper bound for = 98, dropping gradually to 
the lower bound at = 100, while the optimal preventive control w 2 
and treatment control w 3 both maintain the maximum of 100% for 
= 99 before they drop sharply to the lower bound at = 100. Hence, 
this suggests that minimal effort is required to implement the 
application of the insecticide spray on the bleeding grounds of the 
vector and treatment of the infected human than implementation of 
the preventive control of sleeping under mosquito net at a 3 = 10000. 



Figure 19: The plot represents Optimal control w 1 with a 3 =10000 
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Figure 17: The plot represents Optimal control w 2 with a 2 =10000 
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Figure 18: The plot represents Optimal control w 3 with a 2 = 

10000 
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Figure 20: The plot represents Optimal control w 2 with a 3 =10000 
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Figure 21: The plot represents Optimal control w 3 with a 3 =10000 


Conclusion 

In this research article, an optimal control model of malaria 
disease with standard incidence rate was mathematically formulated 
to study the transmission dynamics of the disease in order to 
analyze the optimum control strategy that would be efficacious to be 
implemented to control the disease at a minimal cost. Three control 
functions were introduced to assess and measure empirically the 
efficacy of the use of sleeping under treated mosquito net to prevent 
direct contact and bite from infected mosquito, the application of 
insecticide spray on the bleeding grounds of the vector and the rate at 
which infected humans are treated at each time of the infection. The 
analysis proved that the optimal control strategies considered have 
an optimum and incomparable results on the reduction of the number 
of exposed and infected humans and mosquitoes as compared to the 
model without control as illustrated in the plot of figures for the 
models with and without controls. The numerical examples depict 
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that despite the standard incidence rate, the proposed strategies are 

effective in the reduction of the number of the exposed and infected 

human and mosquito of the disease. 
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